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Visualizing the kinematics of relativistic wave packets 



: 

o . 
o ■ 

(N ; 
Or 

m : 



> : 
o\ . 
■ 

o ■ 

On i 
O ■ 
■ 

o ■ 

^: 
Oh: 

^— > ' 
G ■ 

a • 

s : 
cr. 




Bernd Thaller 

Institute for Mathematics and Scientific Computation, 
University of Graz 
Hemrichstrasse 36, A-8010 Graz, Austria 
(Dated: February 1, 2008) 

This article investigates some solutions of the time-dependent free Dirac equation. Visualizations 
of these solutions immediately reveal strange phenomena which are caused by the interference of 
positive- and negative-energy waves. The effects discussed here include the Zitterbewegung , the 
opposite direction of momentum and velocity in negative-energy wave packets, and the superluminal 
propagation of the wave packet's local maxima. 
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I. INTRODUCTION 

The numerical simulation of the Schrodinger equation 
and the visualization of its solutions has become an im- 
portant part of quantum-mechanical education on all lev- 
els 0, Hi • But, it is almost unavoidable to stumble 
accross strange phenomena, when one attempts a numer- 
ical solution of the Dirac equation. In this article we 
show visualizations of some of these peculiarities and try 
to give an informal explanation. 

The Dirac equation is the fundamental equation for 
relativistic quantum mechanics. Hence, it belongs to the 
most important equations in modern physics. Among its 
big successes is the very accurate description of the en- 
ergy levels of the hydrogen atom. On the other hand, the 
occurrence of several unexplained paradoxes casts doubts 
on its status and interpretation. 

The unexpected behavior occurs even for innocent 
looking initial conditions. A canonical set of initial condi- 
tions for the time-dependent Schrodinger or Dirac equa- 
tions is given by the set of Gaussian wave packets. They 
describe more or less localized quantum states for which 
the product of the uncertainties in position and momen- 
tum is minimal. On the other hand, the set of Gaussian 
initial conditions covers most cases of practical interest, 
because any wave packet can be approximated by a su- 
perposition of a finite number of Gaussian states. 

The motion of Gaussian wave packets according to the 
one-dimensional free Schrodinger equation shows little 
surprises. The example in Fig. ^ shows a nonrelativis- 
tic Gaussian wave packet with average momentum zero. 
Initially, the wave packet is well localized, but it spreads 
during the time evolution. As the wave packet gets 
smeared out, its height decreases, because its norm as 
a square-integrable function must remain constant. The 
wave packet remains a Gaussian function for all times. 

According to the free Schrodinger equation, the aver- 
age position (x) and the average momentum (p) of the 
wave packet obey the rules of classical mechanics. More- 
over, the spreading of the wave packet is independent of 
the average velocity of the wave packet. This spreading 
of the position distribution would also be observed for 
a cloud of classical particles whose density in position 
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FIG. 1: Spreading of a Gaussian wave packet according to 
the Schrodinger equation. 
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FIG. 2: Time evolution of a Gaussian initial wave packet 
according to the Dirac equation. 



space is a Gaussian function, provided that the momenta 
of the particles also have a Gaussian distribution. 

Fig. [21 shows a numerical solution of the one- 
dimensional free Dirac equation. It has the same Gaus- 
sian initial distribution as the nonrelativistic wave packet 
in Fig. ^ yet its behavior is quite different. The rel- 
ativistic wave packet wiggles back and forth, becomes 
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non-Gaussian for t ^ 0, and soon develops characteristic 
ripples. 

This result is so strange that anybody with some expe- 
rience with quantum mechanics (but not with the Dirac 
equation) would first assume that the numerical method 
is at fault. This is a good example supporting the argu- 
ment that the numerical solution of an equation is rarely 
sufficient to understand a phenomenon. This motivates a 
more careful theoretical analysis in order to understand 
the origin of this strange behavior. 

This article contains several black-and-white images of 
solutions of the Dirac equation. Computer-generated an- 
imations showing the whole time evolution can be found 
on the internet [§j • These movies show both components 
of the solution with a color code for the phase of com- 
plex numbers. Thus they reveal more information than 
the simple black-and-white reproductions of the position 
probability density in this article. The forthcoming book 
[3 will be accompanied by a CD-ROM containing a large 
collection of similar movies, animations, and simulations. 
They illustrate and supplement the theoretical exposition 
in order to provide students with an intuitive understand- 
ing that is hard to achieve by studying mathematical for- 
mulas. See the web site 7] for more details and the ideas 
behind the project Visual Quantum Mechanics. 



II. THE DIRAC EQUATION 

In this article, we discuss the time-dependent free 
Dirac equation in one space dimension. We write it as an 
evolution equation in the familiar "Schrodinger form" 



ifi^ ip(x,t) = H tp(x,t), i>(x,0) = ip (x). (1) 

The free Dirac Hamiltonian Hq is the matrix-differential 
operator 



H Q = c oi p + a 3 mc , 



(2) 



where a\ and a 3 are the famous Pauli matrices, and p = 
—ihd/dx. 

The expression for H can be interpreted as a lineariza- 
tion of the relativistic energy-momentum relation 



E = X(p) = \Jc 2 p 2 + m 2 c 4 . 
The square of the Dirac operator is just given by 

tt2 2 2, 2 4 

H = c p +m c . 



(3) 



(4) 



For numerical computations and for the visualizations 
it is advantageous to use units where h = m = c = 
1. These can be obtained from the Si-units by a simple 
scaling transformation. Hence, in the following we use 
Dirac equation in the form 



Instead of o\ and 03 , we could use any other pair of Pauli 
matrices. This would give a unitarily equivalent formula- 
tion. All images in this article would remain unchanged. 

The phenomena explained here also occur in higher 
dimensions, but the one-dimensional situation is much 
easier to visualize. For the Dirac Hamiltonian in three 
space-dimensions, Pauli matrices are not sufficient; 4x4- 
Dirac matrices are needed instead. The book [8| contains 
more information about the Dirac equation in three di- 
mensions. 



III. DIRAC SPINORS AND THEIR 
INTERPRETATION 

A suitable state-space for the solutions of the Dirac 
equation must consist of vector- valued functions 



ip(x,t) 



ll>2{x,t) 



(6) 



because the operator Ho is a two-by-two matrix. These 
two-component wave functions are usually called Dirac 
spinors. The word "spinor" might appear inappropriate. 
In one dimension, all magnetic fields that could possibly 
affect the spin are pure gauge fields. Hence, the two com- 
ponents of © do not describe the spin. Rather, they re- 
flect the appearance of negative energies. (However, the 
doubling of the spinor components in three dimensions is 
caused by the spin.) 

The free Dirac Hamiltonian in momentum space is the 
matrix 



Mp) = 



cp 



cp 



(7) 



which for each pel has the eigenvalues X(p) and —X(p). 
As the Hamiltonian represents the energy of the quantum 
system, one is led to the conclusion that the possible 
values of the energy of a free particle can be positive or 
negative. 

Correspondingly, the Dirac equation has two types of 
plane-wave solutions, which we denote by u pos and it ne g- 
For each pel, 
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U po= (p; X, t) = —= U pos (p) f^F^fp)* 



2ir ncE 



(8) 



where u pos (p) and w nC g(p) are eigenvectors of the ma- 
trix ho(p) belonging to the eigenvalues A(p) and — A(p), 
respectively. Hence, we have 



Hq mp» (p; x, t) = ±A(p) «p™ (p; x, t) 



(9) 



and therefore u pos (u ncg ) is called a plane wave with pos- 
itive (negative) energy. 

In order to be in line with the formalism of quantum 
mechanics, one usually requires that (for each t) the com- 
ponents of a Dirac-spinor be square-integrable, 



d_ 

di 



<ip(x,t) = (-iai-?-+(T 3 )il)(x,t). (5) / \ipj{x,t)\ 2 dx < 00 (for all t and j = 1, 2). (10) 

V ox ) J -00 
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One should be aware that this mathematical requirement 
is closely related to the interpretation of the solutions. 
Our choice makes perfect sense if we agree with the fol- 
lowing tentative interpretation of relativistic wave func- 
tions: 

Suppose ip is a normalized Dirac spinor, then 

/ (|V>i(M)| 2 + hfe(M)| 2 )^ (11) 

J a 

is the probability to find the particle in the interval 
(0.6) C M. 

Here "normalized" means that the above integral from 
—00 to +00 is equal to 1. Similarly, and in complete 
analogy to nonrelativistic quantum mechanics, 

[\\Mp,t)\ 2 + \Mp,t)\ 2 )dp (12) 

J a 

is interpreted as the probability to find the momentum 
of the particle in the interval from a to b. (The hat 
denotes the Fourier transform with respect to x.) This 
interpretation is consistent with the choice of p = —id/dx 
(the generator of spatial translations) as the momentum 
operator. 

The unitarity of the time-evolution generated by the 
Dirac equation guarantees that the normalization of 
Dirac spinors is time-independent. Hence, if there is one 
particle at the beginning, there will be one particle at 
the end. This remains true for arbitrary external fields. 
Contrary to a common belief, the Dirac equation is not 
able to describe pair production within the framework 
of quantum mechanics. For this one needs the formal- 
ism of quantum field theory where many-particle wave 
functions are built from one-particle solutions (and from 
one-antiparticle solutions, which are related to negative- 
energy solutions). 

Square-integrable wave packets can be obtained from 
the plane waves by superposition. This is the same pro- 
cedure that works for the Schrodinger equation. With 
suitable coefficient functions, any square-integrable solu- 
tion of the Dirac equation can be written in the form 

/oc 
(^posO) u pos (p;x,t) 
-OO 

+ ip ncg (p)u ncg (p;x,t))dp. (13) 

The coefficient functions "0p OS and ipneg can be deter- 
mined from the Fourier transform ip(p, 0) of the initial 
function ij)(x, 0) by a projection onto the positive or neg- 
ative energy subspace: 

kr s (p) = p iz $(p> °) = \ i 1 ± ^|r) $(p. °)- ( 14 ) 

Numerical integration of H13|) provides a method to com- 
pute the free time evolution of an arbitrary initial func- 
tion that can be a useful alternative to a finite differ- 
ence scheme. More about the mathematics of the one- 
dimensional Dirac equation can be found in 
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FIG. 3: Space-time diagram and worldline of the position 
mean value for the wave packet in Fig. |21 

For many reasons, including some of the observations 
in this article, the interpretation given above may be 
critisized. Therefore, one should consider this interpre- 
tation merely as a convenient working hypothesis that 
allows us to apply the well-established mathematical for- 
malism of quantum mechanics. 



IV. EXAMPLES OF RELATIVISTIC 
KINEMATICS 

In this section, we describe and visualize three strange 
phenomena shown by solutions of the free Dirac equation. 
Some explanations will be given in the following sections. 

As a first example, we compute the free time evolution 
of the Gaussian Dirac spinor 

^ 0) = (3L) V4 ex P(-^/16) Q • (15) 

A few snapshots of the solution are shown in Fig. [2J This 
image shows the position probability density \ipi{t, x)\ 2 + 
\%p2{t,x)\ 1 according to the interpretation given above. 
Compared to the simple spreading of the nonrelativis- 
tic Gaussian, this wave packet shows a very complicated 
motion. 

The choice of the numerical constants in (|15|l is moti- 
vated by the following consideration. The Fourier trans- 
form of the wave packet l|15|) is 

$(P,°) = Q 1/4 exp(-4p 2 ) (16) 

It describes a momentum distribution that is well local- 
ized in the interval [—3/4, +3/4]. This corresponds to 
a maximal speed less than 3/5, far below the speed of 
light c = 1. Hence, we can be sure that all observations 
have little to do with relativistic effects due to velocities 
approaching the speed of light. 
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FIG. 4: Time evolution of a Gaussian initial wave packet 
with positive momentum. 



Fig. |3| shows a space-time diagram of this solution. 
Here, the density plot visualizes the position probabil- 
ity density as a function of the space coordinate x and 
the time coordinate t in dimensionless units. The white 
curve shows the world line of the average position, which 
obviously does not obey classical (relativistic) kinemat- 
ics. Instead, the expectation value of the position op- 
erator performs a wiggling motion, commonly known as 
"Zitterbewegung" [9j. Apart from the rapid oscillation, 
the wave packet drifts slowly to the right, although its 
average momentum is zero. Nevertheless, it turns out 
that the momentum distribution (as determined from the 
Fourier transform of the wave packet) is still a conserved 
quantity with average momentum zero. 

The second example, Fig. 0] shows the free time evolu- 
tion of the initial spinor 



<P(x,0) = (^-)* 4 exp(-a: 2 /16-i3a;/4) 



(17) 



This spinor is very similar to (|15fl . In position space, 
it is multiplied by a phase factor exp(— i3x/4). In mo- 
mentum space, this means a translation by 3/4. The 
momentum distribution belonging to (|17J) is therefore a 
Gaussian distribution centered around the average mo- 
mentum 3/4. Moreover, the momentum distribution is 
so narrow that all momenta contributing significantly to 
the wave packet are positive. Nevertheless, the solution 
splits into two parts, and the smaller part moves to the 
left (that is, with a negative velocity). 

The space-time diagram of this solution, shown in 
Fig. again shows the Zitterbewegung of the position's 
mean value. But this time, the oscillation quickly fades 
away. We see that Zitterbewegung is sustained only as 
long as the left-moving part and the right-moving part 
of the wave packet have some overlap in position space. 

The third example is shown in Fig. |SJ It realizes a 
wave packet with positive velocities. The initial wave 
packet was obtained as a superposition of a positive- and 




FIG. 5: Space-time diagram and worldline of the position 
mean value for the wave packet in Fig. |1] 
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FIG. 6: Time evolution of a wave packet with positive veloc- 
ity. 



a negative-energy part: 

lp(x, 0) = tp pos (x) + *Pnc g (x), (18) 

The positive-energy part has positive momentum. In mo- 
mentum space, it is defined as 

4> pos (x) = NP pos cxp(-4(p - 4/5) 2 ) Q . (19) 

The negative-energy part has negative momentum, 

&»g(aO = NP neg exp(-4(p + 4/5) 2 ) Q (20) 

In position space, both parts obviously move in the same 
direction. They interfere with each other, thereby caus- 
ing the ripples in the position distribution. 

Fig. shows a space-time diagram of this motion. 
There is no Zitterbewegung at all; the world line of 
the average position is a straight line. Interestingly, the 
peaks of the wave packet move with superluminal speed. 
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FIG. 7: Space-time diagram and worldline of the position 
mean value for the wave packet in Fig. 

This is strange because, at least in principle, local vari- 
ations of the position probability density are observable. 
Let us now proceed with a theoretical analysis. 

V. PARITY AND DIRECTION OF MOTION 

The average velocity of the wave packet <|15|) shows a 
slow drift to the right, which is clearly visible in Fig. [3] 
It is remarkable that the solution is not symmetric with 
respect to reflections at the origin, although the initial 
condition is. Indeed, the Dirac equation is not invariant 
under the replacement 

i/)(x,t) ^il>(-x,t). (21) 

But the replacement x — > —x in the wave function does 
not describe the physical space reflection. Like any 
Lorentz transformation, the space reflection has a part 
that acts on the components of the wave function. The 
correct way to describe the space reflection in the Hilbert 
space of the Dirac equation is the parity transform 

P : ip{x,t) -> a 3 ip(-x,t). (22) 

Hence, the wave packet ()15f) is not invariant under a par- 
ity transform, Pip is not a scalar multiple of ip. 

An example of a parity invariant solution is provided 
by the Gaussian initial wave packet 

^ 0) = (4^) V4 exp(-z 2 /8) Q , (23) 

which has only an upper component and satisfies 
Ptp(x, t) — i/j(x, t) for all times. Fig.|S]shows a space-time 
diagram of the position probability density for this solu- 
tion. The average velocity is zero, and the wave packet 
shows no Zitterbewegung in the expectation value — 
although the position probability density shows ripples 
similar to the first solution. 



FIG. 8: Space-time diagram and worldline of the position 
mean value for the wave packet l|23|l . 

The average velocity of a wave packet is described by 
the classical velocity operator 

v cl (p) = c 2 P H 1 . (24) 

It corresponds to the definition of the velocity in terms of 
momentum and energy, which is familiar from classical 
relativistic mechanics. Note that this relation between 
velocity and momentum depends, in particular, on the 
sign of the energy. For a wave packet with negative en- 
ergy, a positive momentum p thus corresponds to a neg- 
ative average velocity v c \. This is also reflected by the 
phase velocity of the plane waves (JBJ) ■ The phase velocity 
of the plane waves is 

v ph = sign(E)^-. (25) 
P 

Note that the phase velocity is always faster than the 
velocity of light, and even tends to infinity in the limit 
of small momenta. This does not matter, because no 
information can be transmitted with phase velocity. Ac- 
tually, a plane wave is spread over all of space-time and 
all information that is carried by a plane wave is already 
everywhere. The group velocity of relativistic wave pack- 
ets is always slower than or equal to the velocity of light. 
But the sign of the phase velocity carries over to the sign 
of the group velocity. 

For the wave packet l|15|l . Fig. [5] shows the momen- 
tum distributions of the parts with positive and negative 
energy, that is, the functions \^ P os(p)\ 2 and l^neg^)! 2 , 
according to (|T4*jl . We see that the positive-energy part 
has its momentum distribution slightly shifted towards 
positive momenta, whereas the negative-energy part has 
a negative average momentum. For a positive-energy 
wave packet, a positive average momentum means a pos- 
itive average velocity, as usual. But, for a negative- 
energy wave packet, the negative average momentum cor- 
responds to a positive average velocity. Hence, the whole 
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FIG. 11: Momentum distributions of the positive and nega- 
tive energy parts of the wave packet in Fig. |S] 



FIG. 9: Momentum distributions of the positive and negative 
energy parts of the wave packet in Fig. 




FIG. 10: Momentum distributions of the positive and nega- 
tive energy parts of the wave packet in Fig. 2] 



wave packet has a positive average velocity, which can be 
seen clearly in Fig. [3J 

The opposite direction of velocity and momentum in 
the negative-energy part of a wave packet immediately 
explains the behavior of the wave packet l|17|) in Fig. 01 
For this wave packet, the momentum distributions of 
the parts with positive and negative energy are shown 
in Fig. ^3 We see that both parts consist of positive 
momenta. For the smaller negative-energy part, this cor- 
responds to a negative velocity. Hence, the wave packet 
in position space is a superposition of an approximately 
Gaussian wave packet with positive energy moving to the 
right, and a smaller part with negative energy moving to 
the left (also with a positive average momentum) . 

In case of the third example i|18|) . the negative-energy 
part has negative momenta, as shown in Fig. 1111 Here the 
two parts are of the same size and they have opposite mo- 
menta, they have the same average velocity. Hence, the 
parts with positive and negative energy move together 
in position space. The opposite phase-directions of the 
two parts cause the interference effects. At the average 
momentum (p), the initial state contains the superposi- 
tion of the plane waves u pos ((p); x, 0) + u ncg (—(p);x,0). 
In both components of the spinor, this is proportional 
to a superposition of exp(i(p)x) and exp(— i(p)x). This, 
precisely, causes the sine-like distortions of the Gaussian 
position probability density visible in Fig. |6l 

Hence, the interference ripples have their origin in the 



phase interference of the parts with positive and negative 
energy. As the phase moves with a superluminal average 
velocity i; ph = X((p))/{p), also the interference pattern 
moves with superluminal speed. Because it is just an 
interference pattern, it carries no information from one 
point in space to another. Hence, despite the fact that 
local variations in the position probability density are, in 
principle, observable, the fact that these variations move 
at superluminal speed is no contradiction to the special 
theory of relativity. A classical example for this effect 
is the apparent motion of Moire patterns. If two peri- 
odic patterns with slightly different periodicity lengths 
are slowly shifted with respect to each other, the pattern 
of the superposition (the Moire pattern, which is locally 
observable) can move with almost arbitrary speed. In 
fact the recently discussed phenomenon of superluminal 
tunneling also consists in the motion of an interference 
pattern of plane waves (thereby the transmitted wave 
packet changes its shape, but no information is being 
transmitted with superluminal speed) . A local perturba- 
tion of the interference pattern can only propagate with 
a velocity less than c. 



VI. ANALYSIS OF ZITTERBEWEGUNG 

The Zitterbewegung (that is, the shape of the white 
curve in the space-time diagram) is described by the ex- 
pectation value (in the given initial state) of the operator 
(see, for example, @) 



x(t) 



x + v cl {p)t + Z(t). 



(26) 



Here v c \ is the classical velocity operator (|24|l . and Z(t) 
describes the oscillation, 



Z(t) = (2iH )- 1 (e 



2iH t 



l)(aj x -v c i(jp)). (27) 



Note that ca\ is the instantaneous velocity according to 
the Dirac equation, because the time derivative of the 
position in the Heisenberg picture is just 



— x = i[H ,x] = ca 1 . 
at 



(28) 



The initial wave packets i|15|) and (|15|l are eigenvectors 
of a\, belonging to the eigenvalue +1. Hence the initial 
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velocity of these wave packets is +c, hence in Fig. [3] and 
Fig. [SI the slope of the white curve at (x, t) = (0, 0) is 45 
degrees. 

In momentum space, Z(t) is just a multiplication by 
a matrix-valued function of p. The operator Z(t) an- 
ticommutes with Hq- This means that Z(t) maps a 
state with positive energy onto a state with negative en- 
ergy. Suppose that \ip) is, for example, a positive-energy 
state. Then Z{t)\%p) has negative-energy and is thus or- 
thogonal to \ip). Hence, the scalar product (ip\Z(t)\ip) 
must be zero. The expectation value (ip\Z(t)\ip) can only 
be nonzero for wave packets having both positive- and 
negative-energy parts. 

Wave packets located in different regions of momen- 
tum or position space are orthogonal. The operator Z(t) 
is a multiplication operator in momentum space which 
means that it does not change the localization properties 
of a wave packet in momentum space. Hence, Zitterbe- 
wegung can only be significant, if the function i/) pos (p, 0) 
has a significant overlap with the function ?p ncg (p, 0) in 
momentum space. 

In position space, the operator Z{t) is nonlocal, but 
it turns out, that it does not change the approximate 
localization of a wave packet all too much. Hence, if 
a wave packet ip(x) is (approximately) located in some 
region R of position space, then the function Z{t)ijj(x) is 
approximately located in a neighborhood of that region. 
As a consequence, Zitterbewegung is only significant as 
long as the parts with positive and negative energy are 
close to each other in position space. 

Now, let us consider our examples. Fig. shows that 
the momentum distributions of the positive and negative 
energy parts of the initial wave packet l|15|) have a sig- 
nificant overlap in momentum space. This fact does not 
change with time, because the momentum (and hence 
the momentum distribution of the initial wave packet) is 
a conserved quantity according to the free Dirac equa- 
tion. 

Both parts of the initial wave packet have average mo- 
menta close to zero, corresponding to average (classical) 
velocities close to zero. Hence, the corresponding parts 
will remain close together also in position space. Hence, 
for this wave packet, Zitterbewegung is a sustained phe- 



nomenon. (Actually, the amplitude of the Zitterbewe- 
gung vanishes very slowly, as \t\ — > oo. A mathematical 
argument for the asymptotic decay of Zitterbewegung is 
given in y|.) 

For the wave packet in Fig. the momentum distri- 
butions of the positive and negative energy parts overlap 
completely (see Fig. I1(J[I . Hence, we can indeed observe 
Zitterbewegung, but only as long as the two parts occupy 
approximately the same region in position space. Be- 
cause the two parts have opposite velocities, they quickly 
separate, and the amplitude of the Zitterbewegung de- 
creases rapidly. 

For the wave packet in Fig. the parts with positive 
and negative energy move together in position space, but 
there is no Zitterbewegung, because these parts are lo- 
cated in different regions of momentum space. 



VII. CONCLUSION 

The interference effects occur only for those solutions 
of the Dirac equation that are composed of parts with 
positive and parts with negative energy. The origin for 
the ripples in the position distribution lies in the fact that 
the parts with positive and negative energy have opposite 
phase directions. 

Because any square-integrable spinor-valued function 
(in particular, a Gaussian spinor) is, in general, a super- 
position of positive and negative energies, one is likely to 
run into these phenomena when one computes a numer- 
ical solution of the Dirac equation. 

One may argue that superpositions of positive and neg- 
ative energy states are not physically observable. A wave 
packet with only one sign of energies shows no Zitter- 
bewegung and behaves reasonably. Indeed, in quantum 
field theory, one constructs a Hilbert space of many- 
particle electron states from the positive-energy solutions 
of the one-particle Dirac equation, and the positronic 
states are built from charge-conjugated negative-energy 
solutions. A problem with this approach is that the elec- 
tronic and the positronic states cannot always be sepa- 
rated in the presence of external fields. 
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